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1. Introduction 

There has been a growing interest in chiral perturbation theory (ChPT) for the 
quenched approximation of QCD . The motivation for this is the fact that at present 
the quenched approximation is indispensable for the numerical study of QCD. Since ChPT 
is used in order to analyze and extrapolate numerical data, it is necessary to adapt it to 
the quenched approximation. To this end, we have developed a systematic, lagrangian 
technique that can be used straightforwardly to calculate quenched correlation functions 
within ChPT [Q. Quenched ChPT has been applied to calculate several important quan- 
tities, such as masses and decay constants 

However, aside from these practical results, it turns out that the special role of the 
ry' in the quenched approximation leads to a very singular infrared behavior of quenched 
ChPT, which suggests that quenched QCD does not have a chiral limit [^0. Unlike the 
case in the full theory, the mass of the r]' cannot be taken to infinity so as to leave us 
with an effective action which describes only the pseudo-Goldstone mesons. Instead, a 
double pole term proportional to the singlet part of the r]' mass squared shows up in the 
r]' twopoint function and causes the peculiar infrared behavior. 

In this paper, we wish to extend our investigations to the case of partially quenched 
theories. Partially quenched theories are theories in which not all fermions are quenched; 
only for some of the fermions present in the theory will the determinant in the functional 
integral be replaced by 1. (We assume throughout a bilinear fermion action with only flavor 
diagonal terms.) The same method that we developed for studying ChPT for a completely 
quenched theory can also be applied in this case. Our motivation for considering partially 
quenched theories is threefold: 

First, one may learn more about the peculiar infrared behavior by considering what 
happens when only part of the fermion content of a theory is quenched. In particular, if 
different fermion mass scales are present, one might ask how the infrared behavior depends 
on whether all or only some of the fermions with a common mass are quenched. Also, it 
is interesting to know what happens in the unquenched sector of the theory: is a theory 
with n fermions, out of which k are quenched, the same as an unquenched theory with just 
n — k fermions? In the first three sections of this paper we address these questions. 

Second, partially quenched theories arise naturally in the description of simulations 
in which the valence quark masses are not chosen equal to the sea-quark masses. This 
is a not uncommon numerical technique which, for example, allows one to use Wilson 
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valence quarks and staggered sea-quarks. One would like to have a chiral theory for such 
simulations. 

A third motivation comes from staggered fermions themselves. It is well known that 
lattice QCD with staggered fermions describes QCD with four flavors of quarks in the 
continuum limit. In order to use these fermions for simulations of QCD with only two 
flavors, a trick which has been used is to take the square root of the fermion determinant, 
thereby effectively reducing the number of flavors which appear in virtual quark loops from 
four to two. 

This approach seems justified in weak coupling perturbation theory, but of course the 
question is whether it is really a legitimate technique. Certainly, taking the square root 
is not equivalent to formulating a two flavor theory through a functional integral with a 
local lagrangian. In the continuum limit, however, taking the square root of a four flavor 
fermion determinant (with at least pairwise degenerate quark masses) is exactly equivalent 
to quenching two out of four flavors, and our results about partially quenched ChPT with 
n = 4 and k = 2 should apply. This then allows us to test the idea of taking the square root 
within the context of ChPT: if this trick is legitimate, and does indeed lead to a two flavor 
theory, that should be reflected in ChPT. It means that partially quenched ChPT should 
reproduce the results of unquenched ChPT as long as we allow only unquenched quarks on 
the external lines of correlation functions (including operators which excite bound states 
of unquenched quarks). 

Using results obtained in the first part of the paper, we show that this is indeed the 
case. This is a nontrivial test of the trick of taking the square root and complements 
the argument based on weak coupling perturbation theory, since ChPT addresses a differ- 
ent regime of QCD. As a corollary, we present a simple way in which '>l'su{2) correlation 
functions can be computed in numerical simulations. 

2. Theorems 

In this section we state three theorems about partially quenched gauge theories and 
then give the physical arguments which underlie these results. We leave to sections 3 
and 4 the detailed calculations in quenched chiral perturbation theory |^ which illustrate 
the theorems and various corollaries. For the first two theorems, the physical arguments 
actually constitute proofs, and the calculations of section 4 serve as explicit examples. For 
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the third theorem, however, we give only an argument here only for a special case; the full 
proof will have to wait until section 3. 

We first need to establish our notation. Consider a QCD-like theory with n flavors of 
quarks, Qi. The quark masses rrii (z = 1, . . . n) are completely arbitrary. We then partially 
quench this theory by adding k flavors (0 < /c < n) of pseudoquarks (bosonic quarks), 
Qj, as in ref. 0. Note that the limiting cases of complete quenching or no quenching are 
allowed. The masses of the pseudoquarks are fixed to be equal to the masses of the first 
k real quarks, i.e., rrij [j = 1, . . . /c). In other words the first k quarks are "quenched," 
and the remaining n — k quarks are "unquenched." We call this theory the "5't/(n|/c) 
theory." We do not consider the more general, but probably physically uninteresting, case 
where the masses of the pseudoquarks are arbitrary — in that case there can be virtual 
pseudoquark loops that do not cancel completely against real quark loops. Note that, 
if there are degeneracies between some unquenched and some quenched quarks, i.e., the 
quarks of some mass scale are only partially quenched, one is free to choose which of 
the quarks shall be the considered the "unquenched" ones. In other words, if among a 
degenerate quarks, b {b < a) are quenched, one may arbitrarily choose which b quarks are 
to have indices j, with 1 < j < k, and which b — a quarks are to have indices r, with 
k + l<r<n — k. 

A normal, completely unquenched theory with n quarks will be denoted as an ^^SU{n) 
theory;" it is obviously the same as the SU{n\0) theory. 

The full chiral symmetry of the SU{n\k) theory is the semi-direct productS [SU{n\k) 
SU {n\k) f{\(s)U {1) , where the additional U{1) present in U{n\k)L ® U{n\k)fi is broken by 
the anomaly. 

The special case where all the (i = 1, . . . n) masses are equal is called the "degener- 
ate SU{n\k) theory." In section 4, we also examine another special case where the number 
of quarks and pseudoquarks are even (n — > 2n, k 2k) and the masses just take on two 
values: mi = = = ■ ■ ■ = m2n-i = ^-u, and m2 = = rriQ = ■ ■ ■ = m2n = "n%d- We 
caU this the "doublet SU{2n\2kY theory. 

We define the "super-ry'" of the SU{n\k) theory by the interpolating field 



n 



k 




(2.1) 



^ The product is obviously direct when A; = 0. 



4 



The normalization factor c was taken to be 1/ \/n + k {i. e., 1 / \/6 for n = /c = 3) in ref. , 
but has been left arbitrary here, since different normalizations will often be useful. The 
reader may be confused by the fact the super-ry' does not appear explicitly as the difference 
of the quark r]' and the pseudoquark ry', as it does in the corresponding chiral theory, but 
rather as their sum. The reason is that the mesonic fields of a chiral theory directly 
correspond not to q'j^q, but to tr^qq"^^), since it is the first index of the meson matrix 
E which is the quark index. Taking into account the opposite statistics of quarks and 
pseudoquarks, the relative minus sign between quark and pseudoquark r]' would reappear 
in eq. when written in the latter form. 

We can now state the three theorems about partially quenched theories: 
I. In the subsector where all valence quarks are unquenched (z. e., where all valence quarks 
are of type r, where /c + 1 < r < n), the SU{n\k) theory is completely equivalent to a 
normal, completely unquenched SU{n — k) theory. 
II. The super- ry' (with normalization c = — k in eq. is equivalent to the 

rj' constructed in the unquenched sector of the SU{n\k) theory, and is therefore, by 
I, equivalent to the SU{n — k) r]' . "Equivalent" here means that Green's functions 
constructed from an arbitrary number of super- 77' fields and unquenched quarks, will 
be equal to the corresponding Green's functions with the super-ry' replaced by the 
rj' of the unquenched sector of the SU{n\k) theory (or, what is the same, by the 
SU{n—k) rj'). Green's functions which involve the super-ry' and arbitrary combinations 
of quenched quarks or pseudoquarks are not allowed — there is nothing for them to 
correspond to in the SU{n — k) theory. 
III. Quenched infrared divergences ^-H, coming from a double pole in the 77' propagator 
and associated with some quark mass of mass m-, , will arise if and only if the scale rrij 
is fully quenched, i.e., if there is a pseudoquark of mass rrij for every quark of mass 
rrij. In other words, these unphysical divergences arise if and only if 1 < j < /c (so 
that this quark is quenched) and rrij 7^ m^., for all r with A; + 1 < r < n (so that there 
are no unquenched quarks of the same mass) . 

Theorem I is easily established by a simple argument. Since by supposition all the 
valence quarks are unquenched, the only way the amplitudes could "know" about the 
quenched quarks and pseudoquarks is through virtual loops. But the pseudoquarks have 
been chosen to cancel the quenched quarks exactly in virtual loops, so only the unquenched 
quarks can appear anywhere in a diagram. 
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[Instructions for do-it-yourself graphics: 

For Fig. 1(a), draw an oval (racetrack) with its long axis horizontal. Put an 
"x" on each short end (i.e., the extreme right and left edges). 

For Fig. lb), draw an two ovals as in (a) and put them end to end. Put an 
"x" on the left edge of of the left oval and on the right edge of the right 
oval. 

Good going, you did it!] 

(a) (b) 
Figure 1 

Quark flow diagrams for the rj' propagator; (a) is the straight-through" diagram, and (h) 
is the ^'two-hairpin" diagram. Arbitrary numbers of gluon corrections and virtual quark 
loops (if the theory is not fully quenched) are implicit. 

Theorem II also relies on the cancellation between quenched quarks and the pseudo- 
quarks, but this time in valence lines. There are two quark flow diagrams that contribute 
to the Tj' propagator: the "straight-through" diagram (Fig. la) and the "two-hairpin" di- 
agram (Fig. lb). Note that we only specify the valence quark lines in these diagrams; in 
general (for n ^ k) there will be additional virtual quark loops. By the deflnition of the 
super-ry' and the opposite statistics of quarks and pseudoquarks, the pseudoquarks will 
cancel the quenched quarks both in the straight-through and the two-hairpin diagrams 
(in the latter case, the cancellation takes place separately in each hairpin). Only the un- 
quenched quarks survive in each diagram, and they of course will also be the only survivors 
in any virtual quark loops. Thus the contraction of any two super-ry' flelds is the same 
as the contraction of two SU{n — k) rj' flelds. (The choice c = l/yn — k reproduces the 
canonical normalization.) Similarly, the contraction of a super-ry' with some combination 
of unquenched quark flelds is the same as the contraction of an SU{n — k) rj' with the same 
combination, since only the unquenched quarks in the super-ry' will contribute. 

Theorem III relies on the fact that the infrared divergences that have been found 
[|^-|^ in the quenched theory arise from the two-hairpin diagram (Fig. lb), which gives a 
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double pole in the rj' propagator. In section 3, we calculate the propagator in the neutral 
meson sector in partially quenched chiral perturbation theory, and show that the offending 
double poles can only arise when a mass scale is fully quenched. Here, we make this 
result plausible by examining the degenerate SU{n\k) theory. It is not difficult to argue 
that, when there is only one quark mass scale, double poles arise only when the theory is 
completely quenched, i.e., when k = n. 

Consider the propagator of the r]' of the degenerate SU{n\k) theory (constructed from 
the n quarks only, with no pseudoquarks) . In a normal unquenched theory, the two- hairpin 
contribution to the propagator would include diagrams with arbitrary numbers of virtual 
quark bubbles between the two hairpins. The set of all these diagrams, together with the 
straight-through diagram, is a geometric series which can be summed to a simple pole, 
with an r]' mass shifted away from the common meson multiplet mass by the usual singlet 
contribution. In a partially quenched theory, the sum over bubbles is still present, and 
differs only by an overall normalization (coming from the counting of the valence loops in 
the hairpin) relative to the case where only the unquenched quarks are allowed. Therefore, 
the two-hairpin diagram has an "incorrect" normalization relative to the straight-through 
diagram, so the full propagator will not just be a simple pole with a singlet mass term 
added. We can correct for this mismatch by adding and subtracting the proper amount of 
straight-through diagram. The complete ry' propagator will then be a sum of two simple 
poles, one with mass equal to the common meson multiplet mass (from the subtracted piece 
of the straight-through diagram) and one with a shifted mass which includes the singlet 
contribution. Since there are no double poles, there will be no unusual infrared divergences. 
The only exception occurs for k = n, when there are no bubbles to sum. Then the complete 
r]' propagator is just the sum of two terms: a double pole from two-hairpin diagram and a 
single pole from the straight-through diagram 0,^. 

An example which illustrates several of the ideas discussed above is the construction, 
in a degenerate SU{n\k) theory, of an SU{n — k) rj' out of the diagrams for the SU{n) rj' . 
Let both of these particles be described by canonically normalized fields: 




n — k 




(2.2) 
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The straight-through diagram of Fig. la is then clearly identical for the rj'suin) 
the v'su{n-k)^ since the factor of n or n — A; from flavor counting in the loop cancels against 
the field normalization factors. The two-hairpin diagram (Fig. lb) is however normalized 
differently for the v'suin) ''^'su{n-k)i since there are now two flavor loops. To 

get the v'su{n-k) two-hairpin from the v'suin) two-hairpin, one must multiply the latter 
by (n — k)/n. This holds irrespective of the number of virtual loops in the two-hairpin 
diagram. Note that each virtual bubble is normalized the same in both cases, since we are 
always working in a SU {n\k) theory with a net total oin — k flavors in virtual loops. The 
difference in normalization arises only from the valence (hairpin) loops. 

We may thus obtain the correct v'su{n-k) Propagator in an SU{n\k) theory from the 
v'su(n) propagator in that same theory by making a simple readjustment of the relative 
weights of the diagrams. This will prove useful when we discuss staggered fermions. 



3. Proof of Theorem III 

We begin by writing down the lagrangian for the SU{n\k) theory. Define the (n + 
/c) X (n + k) hermitian field $ by 

(3.1) 



where is the nxn matrix of ordinary mesons made from the n ordinary quarks and their 
antiquarks, is the corresponding k x k matrix for pseudoquark mesons, and x is a A; x n 
matrix of mesons made from a pseudoquark and an ordinary antiquark. The unitary field 
S is then defined as 

E = exp(2z$//) , (3.2) 

with / the tree-level pion decay constant. The (n + k) x (n + k) quark mass matrix is given 
by 

Mij = rriiSij , (3.3) 

where, as discussed in the previous section, the masses for z = 1, . . . , n are arbitrary, 
and we take the pseudoquark masses equal to the first k quark masses: run+j = ruj for 
j = l,...,/c. 

The euclidean lagrangian is then 

£ =Vi($o)str(a^Sa'^St) - V2{^o)str{M^ + M^^) 

+ Vb($0)+"t^5(^0)(aM*0)', ^^'^^ 
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where the functions Vi can be chosen to be real and even by making use of the freedom 
aUowed by field redefinitions 0. We choose $o = str($) which corresponds to c = 1 in 
eq. ( |2.1| ). Since the two-hairpin-like interactions between neutral mesons must have no 
dependence on n or /c at tree level, this choice of c guarantees that the parameters in the 
expansion of the Vi are n- and /c-independent at tree level. 

For the purposes of this section we just need the quadratic terms in ( pl^) . We define 



Vi(O) 

V2iO) = v' ^ ^ 



(3.5) 



Note that V5(0) and Vo"(0) are not the same as in ref. [|1 because of the different choice of 
normalization for $o- 

We are now in the position to prove theorem III. In the case of arbitrary quark mass 
the simplest approach is just to calculate the neutral- meson propagator explicitly in tree 
approximation. We work in the basis of the states Ui, i = 1, . . . ,n + k corresponding to 
uu, dd,ss, . . ., and their pseudoquark counterparts. From eqs.( 3.4|) and (lOI), the neutral 
inverse propagator in momentum space is 

= S,j{p^ + Mf)e, + ^e,e„ (3.6) 

where Mf = Svuii/f'^, we have taken a = (it is easy to reinstate later on by the 
substitution fx'^ ^ /i^ + ap^), and is defined by 

-1-1, for 1 < z < n , . 

—1, for n + l<i<n + k' 



It is straightforward to invert ( p.6| ), either by expanding in powers of or by guessing 
the form of the inverse and fixing the coefficients by GG~^ = 1. We have 

p^ + Mf {p^ + Mf){p^ + M^)F{p^) ' ^ ^ 



where 



+ + (3.9) 

1=1 * r=k+l ' 
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The last equality in (|3.9| ) follows from the fact that the pseudoquark masses have been 
chosen equal to the first k quark masses. 

Theorem III now follows by examination of eq. ( |3.8| ) . We show in the Appendix that 
Fip^) has no double zeros, so no double poles in G can arise from F. Therefore, the 
only way there can be a double pole in G is for Mi = Mj, for some z,j (this is of course 
trivially satisfied for i = j), and Mj ^ Mr, for all j between /c + 1 and n. Since Mr, 
/c + 1 < r < 71, are just the masses of the neutral mesons composed of unquenched quarks, 
the latter condition implies that double poles occur at mass Mj if and only if quarks of 
the corresponding mass are completely quenched. This is just the content of theorem III. 

It is instructive to examine eq. ( ^.8| ) in the degenerate limit {Mi = M for all i). For 
k ^ n, we have 

^ 5,jei - l/{n - k) l/(n - k) 

p^ + M^ p-^ + M^ + {n-k)ijy3 ' ^ ' 

This clearly illustrates a result of section 2: that the degenerate propagator (for /c 7^ n) is 
the sum of two simple poles, one with mass equal to the common meson multiplet mass 
and one with a shifted mass which includes the singlet contribution. Note that for k = n 
one sees immediately from eqs. ( p.8| ) and that there are always double poles in G for 



4. Examples 

We would like to demonstrate the theorems in some explicit examples. First, we will 
consider the case of completely degenerate quark masses, with n normal quarks and k 
pseudoquarks. We have calculated the self-energies of the pion and the super-ry' to one 
loop. The euclidean-space pion self-energy is 



=-^{n-k){p^ + ml) I{ml) + ^^m^ I I{ml) - ^ J,^""/^ 

6 [n — k) \ 1 + ■^a{n — k) 



+ A{n-k) {vi'{Q)p' + 



(4.1) 



-m. 



8v l + la{n-k) 

where ml = 8fm//^, and 

d^q 1 



{2n)^q^+m 



2 ■ 



(4.2) 
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For the purpose of this paper we do not need to specify how we regulate such integrals. 
m^' is the mass of the tj' in the SU{n — k) theory, and is given by 

2 l{n-k)^i^ + ml 
1 + ^a[n — k) 

Theorem I is clearly obeyed by this result: the self-energy is a function oi n — k only, 
and therefore is equal to the pion self-energy computed in a theory with n — k normal 
quarks and no pseudoquarks. 

For k < n the above results only have chiral logarithms of the standard type, arising 
from integrals over single pole propagators. On the other hand, for k = n the result is 

This contributes a term to the pion mass which goes like /U^m^logm^, unlike the usual 
mtlogm?.. It is an example of the "pathological terms" previously seen in quenched 
caLlationsI Thi. is a special ca^ of thece. Ill: the. is a pathology as ^ which 
comes from a double pole and arises only for /c = n, in which case the mass scale m is fully 
quenched. 

To demonstrate theorem II, we have to compute the self-energy for the super- 77', $o- 



On the external lines, we choose c = Xj \Jn ~k in eq. (2J.), so that for /c = we are just 



calculating the conventionally normalized r\' self energy in an ordinary unquenched theory. 
(Recall, however, that $0 in the potentials ^^($0) is normalized with c = 1.) For the 
degenerate case, the result is 



S*o(p) = 1 {{n - kf - 1 I{rnl) + I 

n — k \^ 1 + ^a{n — k) 

+ Z!|:M(„ _ U)ml {[{n - kf - 1] n^i) + 

2v \ 1-1- ^a(n — k) 



-Wi\Q){n-k){[{n-kf-l]mll{ml) + {ml,-p^)—^ 

\ 1 -f ^a[n — k) 



a{n — k) 

(4.5) 



^ In other quantities such as (V'V') or fx/ fir such pathologies can lead to actual infrared 
divergences as a quark mass is taken to zero, see refs. [^^. 
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Again, for k < n we see that this is only a function of n — k, and therefore equal to the 
self-energy of the rj' in the SU{n — k) theory. For k = n the normalization c = 1/ \Jn — k is 
clearly inappropriate, and we should multiply ( [4 .51) through by (n — k). The expression then 
vanishes when n = k (note that in that case rriT^ = rrirj'), consistent with our expectation 
that the $o propagator (with finite c) vanishes in the fully quenched theory 0. 

As a final example, we consider the "doublet SU {2n\2k)" theory, in which we have n 
quarks and k pseudoquarks with mass m^, and n quarks and k pseudoquarks with mass 
nid- We will only present that part of the one loop pion self-energy which comes from the 
4- meson vertex proportional to Vi(0), since the full expression is quite cumbersome. All 
our conclusions hold separately for the contribution from each vertex to the self-energy, 
since the parameters multiplying these vertices are free. The result is (here we set a = 
for simplicity) 



" 3/2 J (27r)4^^ 



1 



1 



+ 



(n — k) 



1 



1 



+ 



1 



+ 



+ 



n — k yq"^ + mfj q"^ + m'jy 
2 / o . 1 



n — k 



cos^e 



q^ -|- 



+ sin^6'- 



1 



q^ -|- m_ 



(4.6) 



(The quartic divergence present in eq. ( [4.6|) is cancelled in the total self-energy by a term 
coming from the measure of the path integral.) In this expression. 



4(m„ +md)v _ 
P ~ 2 
1 



{mfj + m%), 



m +-{n — k)iJ. ) ^ \l -{n — k)'^iJ.^ + -{m 

9 ^ 



m 



2 ^,2 
D) ' 



(4.7) 



and 
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The five different masses which appear in eq. (^]^) correspond to the various meson masses 
which appear in the SU{2{n — k)) theory. In the flavor off-diagonal sector, mu and mo 
correspond to mesons of types UiUj and didj with i ^ j, and fn to du or ud. In the flavor 
diagonal sector, there are 2{n — k — 1) yr'^-like mesons with masses mjj and rriD, and two 
other neutral mesons with masses m± due to the singlet-nonsinglet mixing which occurs 
for rriu 7^ rud- is the mixing angle between these two latter neutral mesons. In the case 
that n — /c = 1, the coefficients of the mfj and mj^ poles in eq. ( [4.6| ) vanish, consistent with 
the meson spectrum of the SU{2) theory. 
For k = n, we have 

which is independent of n since results in a fully quenched theory must depend only on 
the valence quarks. Eq. ( [4.9|) agrees for u ^ s with the result we obtained in computing 
the one loop corrections to the kaon mass (c/. ref. [^). 



5. Application to Staggered Fermions 

In this section we will apply some of the results obtained in the previous sections to 
lattice QCD with staggered fermions. In the scaling region, this theory describes QCD 
with four quark flavors, which can be given nondegenerate masses by using nonlocal mass 
terms [^0. If one would like to consider QCD with two flavors, one can use the so-called 
reduced staggered fermion formalism [^^, which however leads to a complex fermion de- 
terminant 0. Also the reduced staggered fermion action does not possess any continuous 
chiral invariance, unlike "normal" staggered fermions. An alternative is to consider normal 
staggered fermions and define a two-flavor theory by taking the square root of the determi- 
nant |[T^] . This corresponds to quenching two of the four flavors. We therefore expect that 
the low energy meson effective theory will be described by SU {4\2) chiral perturbation 
theory. For this to work, the masses need to be at least pairwise degenerate. 

If only a single site mass term is used, the staggered fermion determinant with degen- 
erate quark masses is positive [^. Since the continuum limit is a degenerate four- flavor 
theory with a determinant which is the fourth power of a one- flavor determinant, one ex- 
pects that taking the positive square root of the staggered fermion determinant leads to 
the desired determinant for the continuum two-flavor theory. 
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If nonlocal mass terms are used, the determinant is not positive in general. However, 
the continuum determinant for each flavor is (formally) positive, so one might expect that 
with staggered fermions close enough to the continuum limit, no problem arises in taking 
the square root. 

In this section, we will consider the definition of two-flavor meson operators in the 
mass degenerate two-flavor theory obtained from the degenerate four-flavor theory in which 
the square root of the determinant is taken. Theorem I tells us that we can obtain the two- 
flavor unquenched theory in this way, and that no problems are to be expected from taking 
the square root. For nonsinglet mesons no tuning of the operators is required because one 
may use the same operators as in the four-flavor theory. 

In general, however, one will need to tune the staggered hadron operators in order to 
project out the various continuum hadronic states of interest . For example, tuning will 
be necessary in order to have more than one different quark mass within the staggered 
fermion formalism [§|,|T^. In particular, one expects that the definition of an operator for 



the '>l'su{2) four-fiavor theory will require tuning, even with degenerate quark masses. 

What we wish to show here is that nevertheless two ways exist for choosing a mass matrix 
and a meson operator which do not require tuning of the operator in order to define a pure 
^'su{2) four-fiavor theory. The first method consists of applying theorem II, whereas 

the second method makes use of a peculiarity of nonlocal staggered fermion mass terms. 

We will start by reviewing some facts about renormalization for staggered fermions. 
A general mass term for staggered fermions is given by |^ 

5'mass = mx{x)xix) + m^x{x)E^{x, y)x{y) - m^^x{x)E^{x, y)E^{y, z)xiz) 

X xy xyz 

^l'^f^aPiXiw)Ec,{w,x)Ep{x,y)E^{y,z)xiz) 

wxyz 

~ A 5Z '^^^c.f}-/SXi^)Eaiv,w)E,3{w,x)E-y{x,y)Es{y,z)x{z). 

vwxyz 

(5.1) 

The operator E is defined through 

^mi,x{x)E^{x,y)x{y) = l^mi,Cii{z)[x{z)U^{z)x{z + fi) + x{z + P')Ul{z)x{z)], (5.2) 

xy 2/i 
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where the are certain site-dependent sign factors (c/. ref. 0). m^i, is taken to be 
antisymmetric. 

This mass term leads to the foUowing mass matrix M for the four flavors that emerge 
in the continuum limit: 



M = m + m^^^ + \m^^{-ii^i^) + mli^^^^ + m^^s. (5.3) 

The 4x4 ^-matrices form a representation of the Clifford algebra + = 25^^, 
and are identified with SU{A) flavor generators in the continuum limit. We will denote the 
terms in eq. ( |5.3| ) with scalar (S), vector (V), tensor (T), axial vector (A) and pseudoscalar 
(P) respectively. This expression for M can be derived from the fact that each shift in 
the /U-direction of the field accompanied by a multiplication with C^, corresponds to a 
multiplication of the continuum four- flavor Dirac field by the matrix 



Qix{x)x{x + fi)^ (5.4) 

It can be shown that this form of the mass matrix is stable under renormalization, in 
the sense that the coefficients m, m^, . . . will only receive multiplicative renormalizations, 
one for each tensor structure in eq. ( |5.3| ). This was explicitly demonstrated to one loop, 
and supplemented with more general symmetry arguments, in ref. . Due to the presence, 
in the massless theory, of shift symmetries and a continuous chiral symmetry (the so-called 
t/(l)e symmetry) there are no additive counterterms. Note that the mass matrix M needs 
to be diagonalized in order to determine what the mass eigenstates are. 

Let us first consider the simplest possible mass matrix, by choosing only the single 



site mass m in eq. ( |5.1| ) to be nonzero, corresponding to four degenerate flavors. In this 



case, the simplest operator for an rj'gjj^^-^ will be 



r]'{x) oc x{x)e{x)Ci{x)C2{x + IXsix + 1 + 2)C^{x + i + 2 + 3)x(x + i + 2 + 3 + 4) 

... . (5.5) 
+ sum over all permutations on the directions 1, 2, 3 and 4, 

which in the continuum limit corresponds to the operator '075'*/' [H ? where is a continuum 



Dirac field with four flavor components. In eq. (|5.5| ), the lattice gauge fields are implicit. 



In this basis, an rj'gjj^^^ would be created by the continuum operator 

15 



/cont 
VSU(2) 



^1 on^ 

10 



Vo o/ 



(5.6) 



Clearly, in order to construct a staggered operator with this continuum limit, we need 
the operator S to get a nonzero trace because the r]'su{2) A^'Vor matrix in eq. ( |5.6| ) has a 
nonvanishing trace, and V, T, A and P are all traceless. In addition, we need an operator 
of the type V, T, A and P, since the matrix contains two zero eigenvalues. The fact 
that these operators renormalize differently from S leads to the need to tune their relative 
coefficient. We conclude that with a single site mass term no explicit ti'su{2) operator can 
be constructed in the four-flavor staggered theory without tuning. The only way to avoid 
tuning in this case, is to compute the diagrams for the v'sui'i) (l!!))' adjust the 
relative coefficients of the straight-through and the two-hairpin diagrams as discussed at 
the end of section 2. 

Actually, the special properties of the tensor operator make it possible to construct 
an il'su{2) without tuning in a different way. To discuss this, we will choose an explicit 
representation of the ^-matrices: 



(5.7) 



In this case, it is necessary to choose a mass term of the tensor type. For definiteness 
we choose 



Mo 



/m ^ 

-m 

m 

V 



(5.8) 



-m 



which corresponds again to four flavors with a degenerate mass m. The minus signs can 
be removed by a nonanomalous chiral transformation. The v'sui^) "^^^^ ^^is mass matrix 



IS 



Projecting to SU{2), we get for the '>]su(2) 



^1 
0-100 

10 

Vo -l/ 



75 



(5.9) 
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^1 on^ 
0-100 



Vo o/ 



= -i/;(-Z^1^2 - «6^4)75'0- (5-10) 



Unlike the previous case, this 'n'su(2) A^ivor matrix is now traceless, which allows us to 
write it as a sum of two tensor terms. The minus sign which appears in this equation is 
removed by the same chiral tranformation that removes the sign in the mass matrix, eq. 
( p.8| ). The relevant observation here is that this v'su{2) A^^o^ matrix is not traceless with 
respect to the mass matrix, i.e. tr (Mo(— z^i^2 — ^^3^4)) 0. This follows from the fact 
that the mass matrix defines what the flavor symmetries are (in the continuum limit). If 
ifjL and ifjji transform under SU (4) l ® SU (4) a; as 

i^L VLtpL, i>R VriPr, (5.11) 

the condition for invariance is 



vImVr = M. (5.12) 



With a degenerate mass matrix as in eq. (|5.8|) , the symmetry group is SU{4:) (in the 
continuum limit). If a chiral transformation is performed to remove the minus signs in eq. 
(p78|), eq. ( |5.12D and the trace condition take on their usual form. 

As mentioned above, the ti'su{2) (|5-9| ) is now constructed from two tensor op- 

erators rather than one scalar and one of some other type. Since all tensor operators get 
renormalized in the same way, no tuning is needed here. The price, however, is the use of 
a tensor mass term, which would make this approach awkward for standard simulations. 
Also, in the case of a tensor mass, in general the staggered fermion determinant is not pos- 
itive (c/. the introduction to this section). Using the v'sui'i) readjusting the relative 
weight of the diagrams by hand, as explained in section 2, will be preferable in most cases. 



6. Conclusion 

In this paper our investigations of ChPT in the quenched approximation of QCD are 
extended to theories in which only some of the quarks are quenched. 

The results are formulated in three theorems. Two of them state that in the sub- 
sector with unquenched valence quarks the theory is equivalent to an unquenched theory 
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with the number of flavors equal to the number of unquenched quarks in the partiaUy 
quenched theory. The super-ry' of the partiaUy quenched theory is equivalent to the ry' of 
this unquenched theory. 

The third theorem deals with the existence of infrared divergences due to the double 
pole in the quenched i]' twopoint function . Such divergences only arise if a particular 
quark mass scale is completely quenched. They do not show up in correlation functions 
with only partially quenched or unquenched quarks on the external lines. 

Some one-loop calculations serve as explicit examples of these results. Moreover, we 
apply the n = 4, /c = 2 case to staggered fermion QCD with a single site mass term, in 
which the square root of the fermion determinant is taken in order to yield two-flavor QCD. 
Our analysis shows that this technique is valid within ChPT, and that the super-ry' of the 
SU{4\2) theory (or equivalently the SU{4) r]' with by-hand reweighting of diagrams) can 
be used to measure SU{2) rj' correlation functions without any of the flne tuning which is 
often necessary for staggered fermions. 

Finally, we have shown that another two-flavor rj' operator, not based on the super- ry', 
exists for which no flne tuning is needed if one employs staggered fermions with a so-called 
tensor mass. 
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Appendix 

In this appendix we present a prooiH of the lemma that we used in section 3. The 
lemma states that the function 

1=1 

has no double zeros if all the ct^ are real. The proof will be by contradiction. So let us 
assume that f{z) has a double zero aX z = (3. First, f{z) diverges when z is equal to any 
of the ctj, so we can assume that (3 ^ ai for all i. Now define z' = z — (3. Then f{z'), 
which is given by 



i=l ' 

with a[ = ai — (3, now has a double zero at z' = 0, with all 7^ 0. Note that the a[ are 
not necessarily real, but can have a common imaginary part. From now on we will drop 
the primes on z and ai and assume that f{z) has a double zero at z = 0. / can be written 

as 

where Pk{z) is a polynomial in z of degree k: 



k 

1=1 

If f{z) has a double zero at z = 0, the constant term and the coefficient of the linear 
term in P{z) have to vanish, i.e.. 

Prom the first of these equations one concludes that the common imaginary part of the ai 
has to vanish. Therefore the at have to be real (in other words, the original double zero 

This proof is due to C. Bender. 



i=l 



1 - 



1=1 



k 

E- 



(6.4) 
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would have to be on the real axis) . Substituting the first equation into the second, we then 
conclude that 

fc 

which has no solution for real ctj. This completes the proof that the function f{z) has no 
double zeros anywhere in the complex plane. 
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